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Densely time sampled multi-frequency flux measurements of the extreme BL Lac object S5 
0716+714 over the past three years allow us to study its broad-band variability, and the detailed 
underlying physics, with emphasis on the location and size of the emitting regions and the evolu- 
tion with time. We study the characteristics of some prominent mm-Zy-ray flares in the context of 
the shock-in-jet model and investigate the location of the high energy emission region. The rapid 
rise and decay of the radio flares is in agreement with the formation of a shock and its evolution, 
if a geometrical variation is included in addition to intrinsic variations of the source. We find 
evidence for a correlation between flux variations at y-ray and radio frequencies. A two month 
time-delay between 7-ray and radio flares indicates a non-cospatial origin of 7-rays and radio flux 
variations in S5 0716+714. 
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1. Introduction 

Blazars constitute a unique laboratory to probe jet formation and its relation to radio-to- 7-ray 
variability. The current understanding implies that relativistic shocks propagating down the jet 
provide a good description of a variety of obsei^ved phenomena in AGNs. To provide a framework 
for the observed flux variations, we tested the evolution of radio flares in context of the standai^d 
shock-in-jet model [1 ; 2]. A shock induced flare follows a particular trend in the turnover frequency 

- turnover flux density {S^ - v^) diagram. The typical evolution of a flare in the Sm - plane can 
be obtained by inspecting the R (radius of jet)-dependence of the turnover frequency, v,„ and the 
turnover flux density, Sm [see 3 for details]. During the first stage, Compton losses are dominant and 
v„, decreases with increasing radius, R, while S,„ increases. In the second stage, where synchrotron 
losses are the dominating energy loss mechanism, v„, continues to decrease while 5,,, remains 
almost constant. Both Sm and v„i decrease in the final, adiabatic stage. As a consequence, the Sm 

- Vm diagram is a useful tool to explore the dominance of emission mechanisms during various 
phases of evolution of a flare. 

We report here a radio to 7-ray variability study of the BL Lac object S5 0716-1-714. We tested 
the evolution of radio (cm and mm) flares in context of the standard shock-in-jet model following 
the Sm - Vm diagram as discussed above. We also investigate the correlation of 7-ray activity with 
the emission at lower frequencies, focusing on the individual flares observed between August 2008 
and January 2011. 

2. Multi-frequency light curves 

A broadband flux monitoring of S5 0716-1-714 was performed over a time period between April 
2007 to January 2011. The multi-frequency observations comprise GeV monitoring by Fermi/LAI 
and radio monitoring by several ground based telescopes. The details of observations and data 
reduction can be found in [4]. Fig. [I] shows the 7-ray and radio frequency light curves of the source. 
The top of the figure shows the weekly averaged 7-ray light curve integrated over the energy range 
100 MeV to 300 GeV. The radio frequency light curves are shown in the bottom of the figure. The 
source exhibits significant flux variability both at 7-rays and radio frequencies. Apparently, the two 
major radio flares (labeled as "R6" and "R8") are observed after the major 7-ray flares. 

3. Evolution of radio flares in the shock-in-jet scenario 

In order to test the evolution of the two major radio flares in the context of a shock-in-jet 
model, we construct the quasi-simultaneous^ radio spectra over different time bins as shown in 
Fig. § (a) [see 4 for details] using 2.7 to 230 GHz data. The observed radio spectrum is usually the 
superposition of emission from the two components : (i) a steady state (unperturbed region), and 
(ii) a flaring component resulting from the perturbed (shocked) regions of the jet. The quiescent 
spectrum (Fig. ^ (b) (dotted cui^ve)) is approximated using the lowest flux level during the course of 
our observations. The quiescent spectmm is described by a power law F{v) = Cfy(v/GHz)"'' with 

' time sampling A/ = 5 days 
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Figure 1: Top : GeV light curve of S5 0716+714 during the first ~3 years of the FermifLA^ observations 
from 2008 August to 201 1 January at E > 100 MeV. Bottom : Radio frequency light curves of S5 0716+714 
observed over the past ^3 years. For clarity, the light curves at different frequencies are shown with arbitrary 
offsets (indicated by a "Frequency + x Jy" label). The major radio flares are labeled as "R6" and "R8". 



Q = (0.92 ±0.02) Jy and a,j = -(0.06 ±0.01). We subtract the contribution of the steady-state 
emission from the entire spectrum before modeling. 

We fitted the flare component spectrum using a synchrotron self-absorbed model, which can 
be described as [see 3; 6 for details] : 

c _c ^v^"' l-exp(-T„,(v/v„,)") 
\v,„J l-exp(-T„,) 



where T„, ss3/2^yl — |^ — 1 j is the optical depth at the turnover frequency, 5,„ is the turnover 
flux density, v,„ is the turnover frequency and Ot and ao the spectral indices for the optically 
thick and optically thin parts of the spectrum, respectively (5 ~ v"). 

The evolution of both R6 and R8 flares in the Sm - v,„ plane is shown in Fig. |2| (c) - (d). In the 
standard shock-in-jet model, oc y^' where £, depends upon the variation of physical quantities 
i.e., magnetic field (B), Doppler factor (5) and energy of relativistic electrons [see e.g. 1;3 for 
details]. The estimated e, values are given in Table 

We notice that there is a significant difference between the theoretically expected (from [1]) 
and our calculated e values (see Table |l|). Therefore, the rapid rise and decay of S,n w.r.t. 
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Figure 2: The evolution of the radio spectra: (a) 230 GHz hght curve showing different periods over which 
the spectra are constructed, (b) Results of a single component spectral fitting at time bin "4", the dotted line 
corresponds to the quiescent spectrum, the dashed one to the flaring spectrum and the solid line to the total 
spectrum, (c) & (d) The time evolution of S,„ax vs Vmax for the R6 and R8 radio flares (see text for details). 



Table 1: Different states of spectral evolution and their characteristics 



Flare 


Time 


bin 


^Calculated 


^Expected 


b 


Stage 




JD [2454000+] 








s=2.2, a=l-2 




R6 


1096-1178 


1-4 


-7±3 


-2.5 


0.7 


Compton 




1178-1194 


4-5 








-0.07 


Synchrotron 




1194-1221 


5-8 


10±2 


0.7 


2.6 


Adiabatic 


R8 


1283-1303 


13-15 


-0.9±0.1 


-2.5 


0.4 


Compton 




1298-1345 


15-18 


1.8±0.2 


0.7 


-2 


Adiabatic 



5 oc r'\ B oc R " and Niy) oc y-' 



particularly in the case of the R6 flare (see Fig. ^ rule out these simple assumptions of a constant 
Doppler factor (5). Consequently, we consider the evolution of radio flares including dependencies 
of physical parameters a, s and d following (7). Here, a, s and d parametrize the variations of 
5 R'^, B oc R " and N{y) oc along the jet radius. Since it is evident that the £ values do not 
differ much for different choices of a and s [7], we assume for simplicity s f» constant. For the 
two extreme values of a = 1 and 2, we investigate the variations in b. The two different a values 
give similar results for b. The calculated values of b for the different stages of evolution of the 
radio flares are given in Table As a main result, we conclude that the Doppler factor varies 
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Figure 3: Discrete cross-correlation func- 
tion (DCF) of the 7-ray light curve w.rt. the 
230 GHz radio light curve. The solid curve 
is the best fitted Gaussian function to the 
DCF curve binned at 1 1 days. 
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significantly along the jet radius during the evolution of the two radio fiares. 

4. Correlated mm-gamma-ray variability 

We apply the discrete cross-correlation function (DCF) [8] analysis method to investigate a 
possible correlation among flux variations at radio and y-ray frequencies. In Fig. 3, we report 
the DCF analysis results of the weekly averaged y-ray light curve with the 230 GHz radio data. To 
estimate the possible peak DCF value and respective time lag, we fit a Gaussian function to the DCF 
curve with a bin size of 11 days. The Gaussian function has a form: DCF{t) = ax exp[^^^~^P-], 
where a is the peak value of the DCF, b is the time lag at which the DCF peaks and c characterizes 
the width of the Gaussian function. The best-fit function is shown in Fig. ^ and the fit parameters 
are a = 0.94 it 0.30, b = (67 it 3) days and c = (7 it 2) days. The significance of the correlation is 
checked using the linear Pearson correlation method which gives a confidence level >97%. This 
indicates a clear correlation between the /-ray and 230 GHz radio light curves of the source with 
the GeV flare leading the radio flare by (67 it 3) days. Consequently, the flux variations at y-rays 
lead those at radio frequencies ~1 month time period, which suggests a non-cospatial origin of 
radio and y-ray emission in the sense that /-rays are produced closer to the central black hole. 

5. Summary and Conclusions 

The evolution of the two major radio flai'es in the v„, — S,„ plane shows a very steep rise and 
decay over the Compton and adiabatic stages with a slope too steep to be explained from intrinsic 
variations, requiring an additional Doppler factor variation along the jet. For the two flares, we 
notice that 5 changes as R^ '^ during the rise and as R^-^ during the decay of the R6 flai^e. The 
evolution of the R8 flare is governed by 5 /?" during the rising phase and 5 R^^-^ during the 
decay phase of the flai^e. Such a change in 5 can be due to either a viewing angle (6) variation 
or a change of the bulk Lorentz factor (F) or by a combination of both. The change in 5 can be 
easily interpreted as a few degree variation in 6, while it requires a noticeable change of the bulk 
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Lorentz factor. A similar behavior has also been observed in a parsec-scale VLBI kinematic study 
of the source, which showed that the jet components exhibited significantly non-radial motion with 
regard to their position angle and in a direction peipendiculai^ to the major axis of the jet [9]. Con- 
sequently, a coiTclation between the long-term radio flux-density variability and the position angle 
evolution of a jet component, implied a significant geometric contribution to the origin of the long- 
term variability. This can be probably a result of precession at the base of the jet, which leads to 
twisted and/or hehcal structures. More observations and modeling is required to understand the 
physical origin of these phenomena. A formal cross-correlation between flux variations at radio 
and y-ray frequencies suggests that 7-ray are produced closer to the black hole. The agreement of 
shock-induced evolution of radio flares with a clear correlation between radio and 7-rays is a hint 
for the shock-induced origin of 7-ray emission in the source. 
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